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' SUMMARY 

An analytical study is presented of a laminar liquid film flowing 
along a surface immersed in a gas stream. ThA results of the analysis 
indicate that: (l) The equilibrium film thicpless is directly proportional 

to the square root of the liquid deposition rate; (2) the equilibrium 
film thickness is directly proportional to the square root of the liquid 
viscosity; and (3) the equilibrium film thickness is inversely proportional 
to the square root of the shear stress at the surface. 

The results are intended to show the parameters affecting the thick- 
ness of this film, but are not intended to be quantitative. A calculation 
made from the analysis does appear reasonable, however, when compared 
with the scant experimental data available. The analysis indicates that 
the equilibrium film thickness will be less in regions of high liquid-film 
temperature (reduced viscosity) and high free-stream momentum. 


INTRODUCTION 

The presence of boric oxide on the surfaces of jet engines utilizing 
boron-containing fuels will probably necessitate a new design of some 
engine components in order to obtain the most efficient operation. If 
the thickness of the deposit on any surface could be calculated as a 
function of engine conditions, this new design could be made with greater 
confidence. Until recently, the two problem areas with regard to deposits 
were the primary combustor and the turbine of the turbojet engine. How- 
ever, recent advances in the design of a primary combustor for pentaborane 
fuel (ref. l) indicate that boric oxide deposits in the combustor can 
probably be greatly reduced. Therefore, the turbine is now the major 
problem area.. 

In the turbine, the flow-area relations are quite critical, and any 
change in these relations due to deposits on the blades should be 
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considered during the turbine design. Either the flow-area relations or 
the operating conditions (temperature, pressure, and pressure ratio) of 
the turbine would have to be changed to obtain maximum efficiency with 
deposits. 

The purpose of the analysis reported herein is to determine the 
parameters affecting the equilibrium thickness of boric oxide on jet- 
engine surfaces, so that the engine designer will have an indication of 
the optimum design conditions. This report does not intend to give a 
comprehensive analysis of the equilibrium thickness of the boric oxide 
film. The results of the analysis are semiquant itative and should be 
used only over the range of conditions included in the assumptions listed 
in the ANALYSIS section. Before the analysis, several experimental tests 
were conducted in which motion pictures were made of^a boric oxide film 
on a flat plate. These tests were to determine the characteristics of 
such a film. The results of these tests and the methods by which they 
aided in the analytical model are also presented. 

The analysis is rest^^ted to a laminar liquid film and assumes that 
the boric oxide in the free stream is in the form of liquid particles 
of near molecular- size. 


EXPERIMENTAL APPARATUS AND PROCEDURE 

The apparatus used in the photographic studies of boric oxide film 
on a flat plate was very similar to that reported in reference 2. Tri- 
methyl borate - methyl alcohol azeotrope was burned in a single tubular 
combustor supplied with air from the central laboratory air supply. At 
the combustor exit a transition section reduced the flow area to a rec- 
tangular opening 3.7 by 4 inches wide. The jet emerging from this opening 
was allowed to enter the room and pass over a flat plate positioned on 
the centerline of the jet with the leading edge about l/2 inch from the 
transition section exit. About 8 inches of free space was allowed for 
photographs. The free jet was then captured by a duct leading to water 
sprays and then to the building altitude exhaust system, which was main- 
tained at as low a pressure as possible in order to ensure the exhausting 
of most of the combustion products. A diagram and a photograph of the 
apparatus are shown in figures 1 and 2, respectively. The fuel system 
was the same as that reported in reference 2. 

The instrumentation used in the tests is summarized in the following 
table: 
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Combustion airflow 
Fuel flow 

Combustor -inlet temperature 
Combustor-outlet temperature 
Flat-plate surface temperatures . 

Flat-plate static pressure 


ASME orifice 
Rotating-vane flowmeter 
Average of 3 thermocouples 
Average of 6 thermocouples 
14 Thermo couples, 1 in. apart on 
plate centerline 
Assumed ambient 


A preliminary test was conducted to determine the persistence of the 
potential core of the free jet. In this test, a movable total-pressure 
probe was used to survey the core of the jet. The test was made at 
ambient temperature with a jet velocity of about 800 feet per second. 

The results of this test showed that the total pressure, 1/2 inch above 
the surface of the flat plate, was constant for the entire plate length. 


EXPERIMENTAL OBSERVATIONS 



RESULTS 


Experimental testing covered the following conditions : 


Core temperature, °F 1000 to 1600 

Core velocity, ft/sec 600 to 1200 

Static pressure (ambient), atm 1 


The motion pictures showed that there was no observable change in 
the thickness or characteristics of the plate deposits after 20 to 30 
minutes of operation. In other words, the equilibrium condition had 
been very closely approached. 

The motion pictures and inspection of the flat plate after shutdown 
showed that the liquid film was smooth (referred to as the laminar region) 
over the first 2 to 5 inches (depending on the test conditions), and 
then a transition to a wavy film (turbulent region) occurred. These 
regions are illustrated in the following sketch: 

»*Flow 



Flat plate 
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It was also observed that when a large liquid particle (torn from the 
combustor liner) impinged on the leading edge of the flat plate, the 
film behind the point of impingement was always turbulent. 

Measurements of the laminar film thickness just upstream of the 
transition point showed that it varied from 0.010 to 0.030 inch depending 
on the conditions. The film thickness varied in a manner similar to that 
predicted by the following analysis; that is, the film was thinner at 
high free-stream velocities and temperatures. 

Several determinations were made of the turbulent-film velocity 
about 6 inches from the leading edge of the flat plate. These values 
were taken from the motion pictures and varied from about 0.0004 (gas 
temperature , 1450° F; gas velocity, 650 ft/sec; plate temperature , 

1370° F) to 0.008 foot per second (gas temperature, 1480° F; gas velocity, 
1150 ft/sec; plate temperature, 1380° F). A photograph of the flat plate 
after a typical run is shiCra in figure 3. 

c 


ANALYSIS 

A completely theoretical analysis of the turbulent region of the 
liquid film would be quite difficult if not impossible. However, an 
empirical analysis can be made after experimental data become available. 
Consequently, only the laminar liquid film will be treated analytically 
in this report . 

From observations of liquid film on the flat plate and deposits on 
the stator blades of a two-stage turbine engine operated on pentaborane 
(ref. 3), it was believed that an analysis of the 1 ami nar film would be 
helpful to the designer because in a number of cases it appears that a 
laminar film would be present. An analysis of the laminar film follows. 

It should be again emphasized that the analysis is only intended 
to show the relation among the variables affecting the equilibrium thick- 
ness and is not intended to be a quantitative calculation method from 
which the equilibrium thickness can be obtained. 

The model assumed has a laminar liquid film on a surface immersed 
in a flowing gas stream. Liquid particles suspended in the gas stream 
deposit onto the surface at a known rate. The shear stress caused by 
the. gas stream acting on the liquid film induces a film motion in the 
direction of the gas flow. A point is reached where the volume flow rate 
past any station in the film is equal to the total volume of liquid being 
deposited upstream of that station. This is the condition of . equilibrium. 
Once equilibrium has been reached, provided that shear stress and the 
deposition rate are both held constant, the film thickness will also 
remain constant. 
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For simplification of the analysis, the following assumptions are 

made : 


(1) The rate at which the liquid deposits on the surface is known. 

(2) The film is l a m i nar (free from waves). 

(3) The film moves at a negligible velocity when compared with the 
velocity of the gas stream. Therefore, the shear stress' acting on the 
film will be the same as the shear stress acting on the surface if no 
film was present (verified by the experimental observations). 

( 4 ) The shear stress across the film is constant (in the y- 
direction) . 

(5) The film is at an equilibrium condition. 

Two additional conditions will also be cdAsidered assuming that (1) 
the liquid film is isothermal in the y-direction, and (2) a known temper- 
ature gradient exists across the liquid film in the y-direction. 


Isothermal Condition 

For the isothermal condition it is assumed that the film is iso- 
thermal in the direction normal to the surface. 

Consider an element of the liquid film (abed of the following sketch) 
x-distance from the leading edge of a surface. 


Direction of flow 
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The total volume flow rate of liquid through the element V x can he 
expressed by 

V x = J j Vdx dz 

where V is the local deposition rate. Using the equation of flow 
continuity for the film results j.n 


w x _ 

= = ^f^x^eq^ 2 (2) 

From the definition of shear stress in a fluid (ref. 4), 

or from the use of assumption 4 and assuming no temperature gradient 
(t and are constant), 

/du\ g'tf 

teJ f = nr = constant (3b) 

Therefore, the velocity gradient through the film is linear and 
varies from u^ at the liquid-gas interface to 0 at the surface so that 

% = \ "i (4) 

When equation (4) is substituted into equation (2), 

V x= 2 u i,x" b eq,x z (5) 

The following expression results from setting up equation (3a) for 
integration: 



( 6 ) 


which gives 
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Equations (5) and (7) can be solved simultaneously for either the 
equilibrium thickness t e q or the interfacial velocity ui to give the 
following results : 



Nonisothermal Condition 

Where a large temperature gradient exists across the liquid film 
(in the y-direction) , the same general analysis shown previously can be 
used (with one exception). This exception is (^cussed as follows: 

Equation (3a) can be written 


~ Tfg J 0 ? (io) 

In order to solve equation (10), the viscosity 1 1 must be expressed 
as a function of temperature, and the temperature inserted as a function 
y [T — f(y)> [L — f(T)]. (A curve of viscosity against temperature 
for boric oxide (from ref. 5) is presented in fig. 4.)' The result of 
this substitution is an expression for the local velocity in the liquid ’ 
film u y as a function of y. With this information, the average film 
velocity uf as a function of t e q may be found from 



Equations (10) and (ll) together with the necessary temperature- 
viscosity relations can be used simultaneously with equation (2) to 
■'btain a value of the equilibrium thickness. 
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Application of Equations (8) and (9) to a Flat Plate Vith 
a Turbulent Gas Boundary Layer 


The equation for the shear stress on a flat plate with a turbulent 
gas boundary layer is given in reference 4 as 


_ 7 PO tt 2 d5 
" 72 g U 0 S 


( 12 ) 


where 


dS 0.296 


(13) 


and therefore. 



0.0288 


Po u o 




g 


(14) 


where 


Ren - 


XU QPQ 

PO 


From assumption (3), 


T f = T s 


(15) 


(16) 


Substitution of equations (14), (15), and (l6) into equations (8) 
and (9) results in the following expressions for t e q and ui: 


°eq 


69.5 V xl i z x0. 

“ \, rt 0 . 8 TT 1.8 0.2 

" zp Q ^J 0 PO 


(17) 


and 


"i = 


4 


n nano ^0.8tt1-8 0.2 

G .0576 V^Q °u 0 ^L 0 

0 ? 


( 18 ) 


Now 


v x = 


(19) 
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where r x (which is a function of temperature, pressure, and fuel-air 
ratio only, see ref. 2) is the average deposition rate (ft/sec) over the 
distance x. Substitution of equation (19) into equations ( 17 ) «.n^ (18) 
gives 


V 69.5 fy^ix 

pg-X - 8 ^' 2 


( 20 ) 


and 


u,- = 




0576 r x PQ ,0 Uo' 


FT ). 2 0.8 

M’O *0 


(21) 


which express the equilibrium: thickness of boric oxide and the inter- 
facial velocity in. terms of known variables (r^ can be obtained by the 
method of ref. 2) . 


Equations (20) and (2l) were solved for the condition of an isothermal 
liquid film over a temperature range from 1000° to 2000° F and for the 
following fixed conditions: 


x, ft .' 0.08 

U 0 , ft/ sec 1000 

p 0 , ft/cu ft 0.024 

r x , ft/sec . . 5xl0 -4 

The results of the calculation are presented in figures 5 and 6 . 

The solution indicates that qualitative agreement was obtained (well 
within one order of magnitude) between the analysis and the experimental 
observations discussed previously. No quantitative experimental data are 
available to compare with the theoretical values of figures 5 and 6. The 
experimental values of film thickness obtained previously are merely 
average values for the combined laminar and wavy portions of the film. 

As previously noted, however, the qualitative effects of the operating 
variables on liquid-film thickness appear to be predicted by the theo- 
retical equation. 

Point of Transition from a Laminar to a Turbulent Liquid Film 

Reference 6 shows that for a liquid film, flowing along the walls 
of a tube, transition from a laminar to a turbulent motion was observed 
when the value of y + was between 12 and 21, where 
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y+ = 



( 22 ) 


Equation (22) was put in terms of variables consistent with the 
liquid-film analysis reported herein and was used to indicate where flow 
transition might be expected to occur. 


Substitution of equation (7) into equation (22) gives 


y+ = 


, n/n / 


y 


(23) 


Substitution of t e q for y and inclusion of all the terms on the 
right side of equation (23) under the radical yield 




y + = 


4 


Wl*e 


eq 




eq 


(24a) 


which reduces to 



(24b) 


A film Reynolds number was defined as 


Ref 


teqttjP l 

H 


(25) 


Therefore , 


. Bef = (y + ) 2 (26) 

The following can then be concluded: for laminar film, Re^ < 144; 

for the transition region, 144 < Ref < 441; and for turbulent film, 

Re f > 441. 

For the conditions at which the sample calculation was made (in the 
previous section), the film Reynolds number varies from 1.685*10“5 (at 
1000° F) to 233X10~5 (at 2000° F) . It was observed in the motion 
pictures that transition occurred at a Reynolds number several orders 
of magnitude less than the transition Reynolds number reported in 
reference 6. 
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Discussion and Application of Results 

The purpose of the analysis reported herein was only to determine 
the relation "between the variables influencing the equilibrium thickness. 
It was not the purpose of the analysis to derive an exact quantitative 
expression relating equilibrium thickness with any set of variables. On 
this basis, two general conclusions can be drawn from equation ( 20) : 

(1) The equilibrium thickness is approximately proportional to the 
square root of the liquid viscosity. 

(2) The equilibrium thickness is approximately inversely proportional 
to the square root of the free-stream momentum (PqUq) . 

From these conclusions, the two following design practices can be derived. 
On surfaces where excessive build-up of boric oxide can cause losses in 
engine performance, it is desirable to ( l) keep the surface temperature 
as high as possible (thereby lowering the viscosity); and (2) keep the 
momentum of the free stream (PqUq) as high as possible. 

It should be pointed out, however, that, although a quantitative 
analysis was not the intended purpose of this report, the results of the 
sample calculation appear quite reasonable. 

The exact point at which transition from a laminar to a turbulent 
liquid film occurs is unknown. On the basis of the results reported in 
reference 6, it appears that a laminar film would be present over a large 
portion of the conditions of interest. 

The experimental observation reported herein indicates that transi- 
tion occurs much sooner than reported in reference 6, although a laminar 
film would still be present over a range of conditions of interest. The 
reason for this disagreement is unknown, but may be because of surface 
roughness or impingement of large droplets on the leading edge of the 
plate. 

The experimental observations also point out that a lamirar film 
would not be present downstream of a point where large liquid droplets 
were allowed to impinge. 


CONCLUDING REMARKS 

From the qualitative experimental observations and the analysis re- 
ported herein the following conclusions may be drawn: 
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1. If a laminar liquid film is present, its thickness will he less 
in regions of high free-stream momentum (p Q u2 where p Q is free-stream 

density and Uq is free-stream velocity) . 

2. If a laminar liquid film is present, its thickness will be less 
in regions of high liquid-film temperature. 

3. The liquid film will be turbulent downstream of points where 
large liquid droplets are allowed to impinge. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, June 14, 1957 
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APPENDIX - SYMBOLS 

g acceleration due to gravity, 32.17 ft/sec^ 
r x average deposition rate over distance x, ft/sec 

Ref Reynolds number based on liquid film 
ReQ Reynolds number based on free stream 

T temperature , °F 

t e( j equilibrium film thickness, ft 

Uo free-stream velocity, ft/sec 

u local velocity, ft/sec 

Uf average film velocity, ft/sec 

Uj_ interfacial velocity, ft/sec 

Uy local film velocity y distance from surface, ft/sec 
V local oxide deposition rate, cu ft/sec 

V x integrated film volumetric flow over distance x, cu ft/sec 

w x integrated film weight flow over distance x, Ib/sec 

x distance parallel to axis of flow, ft 
y distance normal to axis of flow, ft 
y + wall distance parameter defined in equation (22) 
z distance across width of surface, ft 
6 boundary- layer thickness, ft 

p viscosity, lb/ (ft) (sec) 

Pj liquid viscosity, lb/ (ft) (sec) 
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Hq free-stream viscosity, lb/(ft)(sec) 

liquid density, lb/cu ft 
Pq free-stream density, lb/cu ft 

T shear stress, lb/sq ft 

Tj. shear stress in film, lb/sq ft 

T s shear stress acting on surface, lb/sq ft 

Subscripts : 
f film 

x distance from leading edge 
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Figure 2. - Experimental apparatus 
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Figure 3. - Flat plate after typical run (run, 64 min). 
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Figure 4. - Viscosity of boric oxide liquid. 




Equilibrium thickness of boric oxide on a flat plate, thousandths of an in. 
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Temperature , °F 


Figure 5. - Equilibrium thickness against temperature. Distance 
from leading edge, 0.08 foot; free- stream velocity, 1000 feet 
per second; free-stream density, 0.024 pound per cubic foot; 
average deposition rate, 5xi0" 4 foot per second. 
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Figure 6. - Interfa 


